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orsional vibrations, due to nonlinear coupling of the 
torsional and lateral modes, are responsible for a 
“snapping” action of the rotor, which can intensify 
crack propagation. These vibrations can also initiate a rotor 
crack in the area of stress concentration, and can be the cause 
of coupling failure. The torsional vibrations are due to non- 
linear coupling of the torsional and lateral modes, excited by 
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unbalance and the radial sideload force together with the 
rotating anisotropy of the rotor. High torsional response has 
the frequency of 1X and 2X, and the snapping action occurs 
when the rotative speed coincides with the lowest torsional 
natural frequency, or half of it (the 1X and 2X torsional reso- 
nances). Characteristics of the synchronous dynamic stiffness 
around the 1X lateral resonance can be used for early detec- 
tion of a lateral crack on any rotor. This article describes a 


study where the torsional and lateral vibrational responses of a 


rote with rotating shaft anisotropy due to a transverse crack 
were modeled. 


Notation 


a M a ee Shalt lateral and torsional damp- 

ing, respectively 
Inboard disk polar moment of 
inerua 
Shalt lateral sciffnesses in the 
“weak” and “strong” directions 
Median stittness 

- Torsional stiffnesses of coupling 
and shalt, respectively 
Mass of the outboard disk and 
racial sideload applied to it 
Unbalance radius 
Stiffness anisotropy factor: effect 
of shalt crack 
Time 
Outboard disk angular (torsion- 
al mode) coordinate 
Outboard disk torsional variable 
Lateral mode natural frequency 
Torsional mode natural 
frequency 

- Inboard disk angular (tarsional’ 
coordinate 
Natural frequencies of torsional 
modes 
Rotative speed 
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Turbomachinery design changes 


The trend to increased efficiency and higher speeds in 
turbomachinery has resulted in the design of more Hexible 
shafts, which operare at rotative speeds above several of their 
natural frequencies. Their new design makes them more 
prone to cracks due to high-cyele and low-cycle fatigue. The 
accompanying increase in the operating temperatures leads to 
the shaft's vulnerability to creep and stress corrosion cracking. 
At least 28 rotor failures due to cracks were documented 
during the previous ten years in lhe US power industry [1]. 
This explains the continuing interest of scientists and engi- 
neers in the dynamic behavior and early detection of cracked 
rotors, 


Many published papers have discussed what is considered 
the major effect ol cracked rotors, the lateral response [1-12]. 
By comparison, the number of papers on torsional vibration 
due to cracks is significantly smaller [13-15]. Yet, torsional 
vibrations, especially in resonance conditions, are more sensi- 
tive to system dynamic changes because the torsional mode 
damping, which comes from the shaft material internal fric- 
tion, is very low, 


Of the numerous early cracked shaft cases we studied while 


monitoring rotating machines in the field, more than 70% of 


the carly warnings have been in synchronous vibration (1A) 
changes. There is a simple rule, ‘lf the shaft cracks, it is usually 
bowed, and thus produces a new unbalance situation.’ The 
remaining 30% showed activity at the twice rotative speed 


motor 


= 





(2X). Shaft-observing XF probes are the best transducer 
method to monitor rotors. For both 1X and 2X crack symp- 
toms, ether phase or amplitude changes (not always increas- 
ing amplitude} are important. The concept of the 
‘snapping (accelerated) action when the rotor passes over 
top dead center, pioneered by the first author of this paper, as 
q major lateral action, results in the primary coupling from 
torsional-to-lateral vibrations. The concept of “breathing and 
gaping cracks,” while fine in theory, docs not seem to appear 
in practice. In field cases of cracking shafts, the crack appears 
to be continuously open. 


This article is a continuation of the authors’ study of rotor 
lateral‘torsicnal cross-coupling, due to the shaft crack, unbal- 
ance, and radial sicleload force. Special attention is paid to the 
2X torsional excitation, which is responsible for the rotor 
“snapping” action, which may intensify crack propagation. A 
computer simulation of the system enabled the authors to 
investigate a range of rotor parameters, and the parameters’ 
influcnce on the rotor lateral and torsional responses. 


Mathematical model 

The system includes the motor, which is modeled as a 
source of the regular rotation with rotative speed a, con- 
nected to che main shaft through the coupling with torsional 
stiffness and damping K. and D.. The main shaft carries two 
disks and is supported in two relatively rigid bearings. The 
inboard disk with the polar moment of inertia, Fa, has only 
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Figure 1 
The rotor with simulated transversal crack. 
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one degree of freedom, described by the angular coordinate, œ 
(phi). The driving torque is applied to the inboard disk from 
the motor through the coupling. 


The inboard disk is connected to the outboard disk through 
the shaft, which has a torsional stiffness, K, and lateral stiff- 
nesses K Ko (K, < K- The. shaft stiffness anisotropy is 
caused by added local fle: xibility due to a crack, The outboard 
disk has a mass, M, ancl polar moment of inertia / = i, + M, 
where r is the distance between the geometric and mass 
centers of the disk. The radial sideload, P, is also applied to 
this disk. The outboard disk motion is described by one 
torsional coordinate, y (gamma), and two lateral displace- 
ments: x is horizontal andy is vertical in the stationary coordi- 
nates. Nore that when y =0 or r , the shaft “strong stiffness” 
axis y passes through the direction of the radial sideload, P. 
This creates the condition for rotor “snapping” (accelerated) 
torsional action with all possible damaging consequences. The 
system has lateral damping, D, 

When all cross coupling terms are neglected in the original 
four-degrec-of-freedom model, the rotor torsional modes 
have the following natural frequencies: 


ye (2K, 4K - KOLAK? ) 
I} 
Two particular cases, which are N for technical 
applications, namely, soft and stiff coupling, are discussed 
below. 
2.1) Saft coupling: KK >I 
In this case, the first torsional mode represents the twist 
angle between the inboard and outboard disks, and has a 
natural frequency of approximately Y 2K, / 7, . The second 
mode represents an angular motion of the main shaft as the 


rigid ooo: and has a natural frequency of approximately 
Y Kf (2i) - This is less important in practical applications. 


2.2) Stiff coupling: K, K << | 

In this case, the first torsional mode represents an angular 
displacement of the inboard disk, and has a natural frequency 
of approximately VK L, The second mode represents an 
angular displacement of the outboard disk, and has a natural 
frequency of approximately VK il, fw w >w a). For fre 
quencies below the first natural frequency , the angular coordi- 
nate of the inboard disk describes the regular rotation 
(=R. The investigation below considers this particular 
Case, 


I 
=— (2K, +K 44 KITAR? 
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Figure 2 
Overall amplitude of the rotor torsional response 
obtained from the computer simulation of case (a). 


Numerical simulation and analysis 


The computer model of the described rotor was built. Three 
sets ot parameters used in the numerical simulation are pre- 
sented in Table 1 (sce Notation): 


since the numerical integration was performed in the range 
of rotative speeds from 0 ta 5000 rpm, case (a) obviously does 
nat experience any significant lateral motion activity, except 
static lateral displacement due to the radial sideload. The only 
important equation for che outboard disk torsional vibration is 
simplificd to the following: 


. DO. K, qe gPe 
i) + abt y= T cosh yey) + isin [2 + they ghi 
Fy fy T Ri y (2) 


where y, is a stationary component of the angle y, 
Y= yh + Vp, 

Fy. (2) can be considered as a linear oscillator with two 
forcing terms: 1X excitation due to the combination of unbal- 
ance, shaft stiffness asymmetry, and radial sicdeload, and 2X 
excitation due to the shaft stiffness asymmetry and radial 
sideload. The latter excitation results in the effect called the 
“snapping” action; the torque experiences two positive peaks 
when the strong stiffness axis passes the radial = direc- 
lion with the angle 45° (Le, y= TYE ¥ g= rj ais AD. The 


results of the numerical simulation for case {a} are Ka 
in the form of overall torsional amplitude versus rotative 
speed (Fig. 2), Note that the effect of the synchronous excita- 
don is relatively small due ta the low amount of unbalance. 1X 
and 2X filtered torsional responses are shown in Figure 3, 
indicating that two peaks in Figure 2 are, indeed, 2X and 1X 


resonantes, 
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Table 1 
The Main Parameters of the Rotor System Numerical Simulation. 
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Bode plots of the torsional response of the rotor for case (a). 
Aj 1X filtered response, B) 2X filtered response. 
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Figure 4 
The waveform of the rotor torsional acceleration 
at {2 = 1000 rpm for case (a). 
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Overall amplitudes of the rotor horizontal and 
vertical responses for case (b). 
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Figure 6 

Bode plots of the horizontal response of the rotor for case (b). 
A) 1X filtered response. B} 2X filtered response. 
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The torsional acceleration timebase waveform (Figure 4) at 
the rotative speed Q = 7000 rp (close to the 2X resonance at 
@= 1250 rpm) illustrates the fact that, due to the snapping 
action, the tarque experiences two positive peaks when 
y= I and = - The Keyphasor® (once-per-turn marker) 


4 
position here corresponds to y= when the strong axis 
coincides with the direction of the radial sideload. Note that, 
since the waveform was taken below the 2X resonance, the 
torsional acceleration lags the exciration by m, and has nega- 
tive peaks while the snapping torque has positive ones. 


Case {b} (Table 1) is, in a sense, an opposite to case (a); in 
the rolatlive speed range (0 to 5000 rpm), due to stiff coupling, 
the rotor does not exhibit any significance torsional activity. The 
rotor lateral response for this case, obtained by the computer 
simulation, is shown in Figure 5 in che form of overall vertical 
and horizontal vibration response amplitudes. Together with 
LX and 2% Bode plots in Figure 6, they confirm that, in the case 
of a significant difference in the lateral ancl torsional natural 
frequency values, there is almost no interaction between the 
modes, and the lateral response consists of mainly 1X and 2X 
components. The Bode plots of 1X and 2X filtered vertical and 
horizontal responses of the rotor are almost identical 
(Figure 5). 

Case (c} overlaps the first two cases. Both lateral and tor- 
sional natural frequencies are in the range of the considered 
rotative speeds; therefore, the case is expected to exhibit both 
lateral and torsional vibrational activities. The results of this 
case simulation are presented in the form of lateral and tor- 
sional overall amplitudes versus rotative speed (Figures 7 
and 4). 

Note that the overall lateral amplitudes for case (c) are 
practically identical to those of case (b). The same holes trug 
for 1X and 2X filtered vertical and horizontal responses. (They 
are not displayed here.) 
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Figure 7 
Overall amplitudes of the rotor horizontal and 
vertical responses for case (c). 
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Figure & 
Overall amplitude of the rotor torsional response for 
case {c}. Compare with Figure 5. 
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As is shown in Figure 8, the situation is now different for the 
rotor torsional response. In addition to the amplitude peaks 


n v | — 
occurring at 9 = = ande, which were present in Figure 2 for 


case (a), the torsional response of the rotor for case (0 exhibits 
1X and 2X lateral resonances at M=w2 and > respec- 


tively. The existence of fractional combination resonances, due 
; r ; : vere F 

Lo nonlinear coupling, is also shown at Q = a and w+ = 
Figure 9 presents Bode plots ofthe 1X and 2X filtered torsional 
responses [or case (c). The peaks at =v ande, are due to the 
1X lateral and torsional resonances, respectively, and overall 

à i P F 
peaks at Q= = and = are due to the 2X lateral and torsional 


=— 


resonances. 





Figure 9 
Bode plots of the torsional response of the rotor for case (c). 
A) 1X filtered response. B) 2X filtered response. 


Analytical consideration 

The analysis of the full system of four lateral/torsional non- 
linearly coupled equations of motion shows the possibility of 
the sequence of resonances, as is summarized in Table 2. 

As is shown in Table 2, the rotor system without stiffness 


Resonance frequency relations 
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asymmetry fg = @) allows for main lateral and integer combi- 
nation lateral and torsional resonances due to the unbalance. 
The stiffness asymmetry, g, by itself compliments torsional 
resonances and fractional combination resonances of the first 
group. Its product with racial sideload affects all the reso- 
nances of the system. According to Figure $, the torsional 
response of the rotor for case (0 shows fractional combina- 
tian resonances of group 2, which are due to the product 
g PIK. 

As was mentioned in the introduction, in 70% of successful 
early detection of cracks, the warnings were received from a 
1X lateral rotor response change. That makes the investigation 
of the stiffness asymmetry influence on the 1X lateral response 
the most important. The mathematical analysis of the lateral 
motion equations in the region around 1X lateral resonance, 
supported by the numerical simulation, provides the results 
presented in Figure 10, They are in the form of a 1X Bode plot 
and components of the system synchronous dynamic stiffness, 


A} m} 
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Figure 10, 
A) Bode plot of the rotor 1X lateral response in the range of 
rotative speeds close to the lateral natural frequency for the cases 
of q =0,0.1,0.2. 
B) Nondimensional direct and quadrature synchronous dynamic 
sliffmesses in the range of rotative speeds close to the rotor 
lateral natural frequency for the cases g = 0,0. 1,0.2. 


Terms responsible for 
| the resonance are 
proportional to: 


Names of the resonances 


Main lateral resonance 


Fractional lateral 
resonance 


Tarsional resonances 


Integer combination 
resonances 


Fractional combination 
resonances (props |) 


Fractional carnkinatvarn 
TESS [grup 2] 


Table 2 
The list of possible lateral and torsional resonances 


10 Orbit 





June 1995 


It is obvious that the synchronous dynamic stiffness 
depends on the unbalance angular orientation related to the 
shaft weak and strong stiffness directions. Note that in Figure 
10 the change of the stiffness asymmetry factor, g, from 0 to 0.2 
has a very small effect on the phase lag, but the amplitude peak 
changes noticeably, 


As follows from Figure 10B, the stiffness asymmetry primar- 
ily affects the quadrature dynamic stiffness in the region close 
ta the main lateral resonance. Since the synchronous dynamic 
stiffness is the reciprocal of the influence vector used in the 
balancing, the information about its behavior is often avail 
able, and can be used in diagnostics of rotating machine 
health, and early detection of rotor cracks in particular. 


Closing Remarks 

A simple rotor system with a shaft transverse crack was 
considered in this article. The lateral and torsional modes 
interact with each other, due to the rotor unbalance, radial 
sideload, and rotating asymmetry. In the absence of any 
tarsional excitation from the driving and load torques, the 
rotor torsional vibrauiens show noticeable IX and 2X resa- 
nances, This effect was confirmed experimentally [15], The 
main source of the torsional excitation was related to the 
torsional snapping action due to the radial sideload. The 
snapping action can be defined as a periodic (twice per rota- 
tion) elastic energy storage occurring when the shalt strong 
axis approaches the radial sideload direction, followed by its 
release when the rotor weak stillness axis passes through the 
radial sideload direction, 





In addition to the above source of 2X excitation, 1X exciting 
torque also exists, due to the combination of the rotor unbal- 
ance, stiffness asymmetry, and radial sideload. Due to coup- 
ling, the torsional response exhibits 1X and 2X lateral 
resonances. The existence of additional fractional combina- 
tion resonances was numerically found and analytically 
explained, 

The torsional response timebase waveform observation in 
the range around 2X torsional resonance (Figure 4) allows for 
an estimation of the crack angular orientation. The torsional 
response also reflects LX and 2X lateral resonances. The lateral 
response of the system was not noticeably affected by the 
torsional vibration, and consisted primarily of IX and 2X com- 
ponents. (The lateral mode carries more damping.) The syn- 
chronous dynamic stiffness components in the proximity of 
the 1X lateral resonance revealed that the shape of the quadra- 
ture component is most sensitive to the shalt stiffness asym- 
metry. This information can be used for shaft crack detection. 
The latter conclusion can easily be applied directly in the fielel, 
since the information on the synchronous dynamic stiffness is 
usually available from influence vectors used in balancing. 

Rotating machinery designers are normally cautious when 
lateral and torsional resonance frequencies aren't close to 
each other and when the operating speed isn't clase to cither 
the lateral or torsional resonance frequencies or their half 
values, The cracked shaft lateral and torsional response 
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described in this article strongly suggests not having lateral 
resonance Frequency close to half of the torsional resonance 
and vice versa. IF one of these conditions is violated, the 
pending crack will lead to high amplitudes of both lateral and 
Lorsional vibrations, as a consequence of the high level of 
normal and shear stresses in the rotor. The torsional reso- 
nance amplitudes are usually much higher than the lateral 
ones, due to very low damping in the rotor torsional mode. 
Any system changes due to a shaft crack, therefore, would be 
better reflected in torsional than lateral vibration changes. W 


Glossary of terms 


Anisotropic - having properties, such as stillness or conductivity, which vary 
according to the direction in which they are measured. 

Combination resonances - ihe coincidence of excitation frequency with 
the integer and fractional combination of two or more natural frequencies of 
I he resto SYA LETTI. 

Torsional bration - time variation of the angle of pwist, ofa shalt or beam, 
typically measured in tenths of degrees pp. 


Transversal crack -a crack which is perpendicular to the longitudinal axis 
cf the rotor. 

1X- In a complex vibration signal, the notation for the signal component that 
has frequency equal to the rotative speed. Also called synchronous. 


2X, 3X, ete. - In a complex vibration signal, the notation for signal compo- 
nents having frequencies equal to exact multiples of shaft rotative speed. Also 
called higher harmonics, superharmonics. 
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